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A technique is presented for the estimation of airplane longitudinal stability 
and control derivatives over a broad range of angle of attack using data from a 
single large aiqp^litude longitudinal maneuver* the application of a modified stepurise 
regression algorithm to both a complete data string and to that saf»e data partitioned 
into bins as a function of angle of attack is demonstrated* Results from the large- 
scale maneuver agree well with results from 20 maneuvers in which small perturbations 
tram trim at various angles of attack were induced, 

lOTRODOCTION 

ihere has been considerable interest recently in the poststall and spin flight 
regimes of airplanes* Ihe inherently nonlinear nature of the aerodynamic forces 
and rnoments of these flight regimes create difficulties for both the design of 
flight test programs and the application of parameter extraction algorithms* though 
interest and attention have been high, no consistently reliable method of 
metric model identification %ihich can be applied to the transient data collected from 
stall/spin/high-a (angle of attack) maneuvers has been developed* It is ing>ortant 
to be able to build an airplane mathematical model that encompasses all operating 
regimes* Hence, one must find a method of extracting information from the short- 
lived stall/spin/high-a maneuvers in order to confirm a given model or build a new 
model for predicting airplane behavior. 

The Identification of airplane stability and control parameters across a broad 
range of flight conditions can be tedious and costly in manpower, con 5 >uter time# and 
in-flight time* Moreover# parameters cannot be identified by traditional methods in 
all flight regimes. In particular, near stall (both prestall and poststall) regltnes 
and spin entry regimes are often too transient (short lived) to provide enough data 
for analysis by existing techniques. This paper offers a new application of the 
stepwise regression technique that promises to be useful in previously nonidentifl- 
able flight regimes. 

The usual method of parameter identification from flight test consists of first 
trimming the test airplane to some given equilibrium fllglit condition (such as 
a « ^ ^ * Then one or more control surfaces (viz, 

elevator, nadder, aileron) are activated slightly but sharply from their trim posi- 
tion and back to that position. The resulting motion of the airplane is recorded 
along with the control rK>vement. The data string for such a maneuver is usually 5 
to 15 sec long. At data rates of 20 points/sec, the experimenter will then have 
100 to 300 points to analyze. 

Unfortunately, It is not possible to trim a given airplane at all angles of 
attack. Moreover, the passage through certain flight regimes is so fleeting that 
only a few data points are available. In regions where the airplane can be trammed 
and the typical maneuvers can be performed# there is a need to know which of 
the regions provide the most interesting and perhaps nonlinear behavior and hence 
require the closest scrutiny. 

These problems can be solved, as this paper demonstrates# by judiciously using 
the data obtained trom several large amplitude maneuvers. It is shown that for large 



amplitude longitudinal maneuvers initiated from a * one can discern informa^ 

tion on the behavior of the airplane in several regions of angle of attack. In this 
way transient behavior can be analyzed as to tentative model structure and parameter 
values. 

Though the linear model of the airplane is sufficient for small perturbations 
from equilibrium! conditions at low angle of attack, there is little kno%m of the 
proper aerodynamic model at higher angle of attack. Klein has shown in reference 1 
that at high lift coefficient the linear model is inadequate for a general avia- 

tion airplane. In the general case of large excursion in the state variables or 
coupling between angle of attack and lateral motion, Thomas has postulated a need for 
a nonlinear aerodynamic model. (See ref. 2.) 

The maneuvers treated in this paper involve, by their intent, large excursions 
in a. Such large excursions then encon^ss regions of both high and low Hence 

part of the identification process required in this work is that of model structure 
determination. 

The determination of model structure consists of two beric requirements: 

1. Is a linear model adequate? 

2. If the linear model is not adequate, which nonlinear terms must be included 

to provide an adequate model? 

Requirements 1 and 2 are addressed by a modified stepwise regression procedure. The 
stepwise regression is used to enter variables into a regression equation one at a 
time. The modified stepwise regression considers only the linear model terms first, 
after which all candidate nonlinear terms are considered. Upon entry of a nonlinear 
term, superfluous linear terms may be deleted from the model. 

The first section of this report deals with the stepwise regression as it is 
applied to the airplane identification problem. Next, the results of the data anal- 
ysis are discussed in detail. And, finally, a concluding section summarizes the data 
analysis and interprets some of the programmatic aspects of the technique used. Tab 
equations of motion are presented in an appendix. 


SYMBOLS 

b span, m 

Cj lift coefficient, L/qS 

rolling-moment coefficient 
pitching-moment coefficient, M^/qSc 

I 

defined in appendix 
yawing-monent coefficient 
Cy^ longitudinal-force coefficient, 

Cy s^ de-force coefficient, F^/qS 
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normal-force coef f icient , F /qS 

Z 

mean aerodynamic chord, m 
expectation value 

partial statistical F-value, standard errors 

force along X, Y, and 2 body axis, respectively, N 

2 

acceleration due to gravity, m/sec 

moment of inertia about X, Y, and Z body axis, respectively, 

2 

product of inertia, kg-m 
error metric 
lift, N 

mean square due to regression 
pitching moment, N-m 
mass , kg 

number of data points 

body axis roll rate, rad/sec 

body axis pitch rate, rad/sec 

dynamic pressure 

body axis yaw rate , rad/sec 

correlation coefficient between x^ and y 

correlation coefficient between x, and y given that x., 

X K 3 

P and Xp are already in regression 

2 

wing area, m 
ith time interval 

velocity along X, Y, and Z body axis, respectively, m/sec 
airplane total velocity, m/sec 

N X P matrix of N measurements of P model variables 
ith independent model variable 


kg^ 


ijth element of X 



y dependent variable 

-► 

y N X 1 vector of N measurements of dependent variable 

a angle of attack, rad or deg, as appropriate 

a mean value of angle of attack for a bin, deg 

a. value of a at ith measurement, rad 

1 

3 angle of sideslip, rad (6^ produces nose--up pitching moment) 

Aa. -a, - a , rad 

X X o 

6 elevator deflection, rad (standard sign convention) 

e 

-► 

£ N X 1 error vector 

9 pitch angle, rad 

? P X *1 parcuneter vector 

9^ ith element of 0 

p air density, kg/m^ 

^ bank angle, rad 


Subscript: 

o trim value 

Supercripts: 

T transposed matrix 

• derivative with respect to time 
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MODEL SELECTION AND PARAMETER ESTIMATION 

A stepwise linear regression was applied to the data for both model selection 
and the estimation of the parameters comprising the model* A general discourse on 
regression (equation error) techniques is found in reference 3. Itie stepwise pro- 
cedure employed for this report is basically that given on pages 178 through 195 in 
reference 4. In general, a probabilistic parameter estimation method such as maximum 
likelihood is preferred to a regression method. However, as part of the analysis for 
this report, the data were partitioned, destroying its continuous dynamic properties 
over more than a few points* Hence, such probabilistic methods were not applicable* 
Moreover, if rhe model be correct and model variable measurement error small, the 
regression technique gives reasonabl 3 results* (See ref* 1*) 

System identification is divided by Zadeh (ref* 5) into several major 
areas* Two areas that this report considers are model selection and parameter 
identification. 

In the standard linear regression technique, a measured quantity, y, is 
expressed as a function of several independent variables, x^ (i = 1, *.*, P), and 
multiplicative coefficients, 6 ^ ( i - 1, * * • , P) , as 


y = 01 + > 63 X 3 + ... 9pXp 


( 1 ) 


Equation ( 1 ) is derived in the appendix where y corresponds to an aerodynamic force 
coefficient or moment coefficient; the variables x^ correspond to the model 
variables ot, q, 6 ^, and their combinations; the coefficients 9^ represent the 
aerodynamic stability and control derivatives such as 

and • If measurements are made on all variables at times t^, t 2 # t^# 

then equation (1) can be written as N independent equations: 
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which can be summarized in the matrix equation: 

= x5 


-► 

y 


( 2 ) 


-► 

where y is an N x i vector of the N measurements of the dependent variable, 

X is an N p matrix whose jth column contains the N measurements of the 
jth model variable, and o is a P ^ 1 parameter vector (unkno%m) • 

Because of error associated with the measurement of y^, the linear combination 
of model variables may not faithfully represent the system* Hence, one must associ- 
ate an N X 1 error vector e with equation (3> to give 


y = 


X? + e 


(4) 


where C is referred to as the^ "equation error" and often lends that name to the 
linear regression method. If e is stationary white, zero mean, and the ele- 
ments of X are measured without error, then the regression technique gives 

unbiased parameter estimates 0 = [X^X] X^y when the cost function 

J = = (y - X?)’’(y - X^) is minimized with respect to That is, 

E(c X) = 0 is desired. To achieve this, one must have the correct model postu- 
lated in equation (1). Usually, for small longitudinal perturbations, a single 
linear model, incorporating only the model variables a, q, and 6^, will suf- 
fice. However, for larger perturbations from equilibrium flight and/or for high 
it has been shown that certain nonlinear combinations of the model variables 
might be required. (Ser ref. 1.) 

Candidate combinations of the model variables are listed in table 1. The candi- 
date model variables consist of the linear terms a, q, and 6^ and certain 
nonlinear terms. The first three nonlinear terms are aq, and a6^ which yield 

parameters corresponding to the slope of the linear parameters with respect to a. 
Next, two terms which account for symmetrical coupling of the longitudinal motion 
with sideslip (viz, and are consid* *ed. Finally, since the longitudinal 

motion is generally most dependent on a, the hxgher order terms ci^, ..., a® are 
included* If all the candidate variables listed are included, another problem is 
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created: the matrix is singular and § does not exist* *Iherefore, one must 

include enough variables to postulate the correct tftodel E(t*^X) * ^ but not so iRany 
variables that [x'^X] is made singular. 

It is desirable then, to add variables to the model, one at a time, until the 
model is correct* A stepwise regression does exactly that* In a stepwise regres^ 
Sion, the equation error is decreased in "steps" by adding the independent model 
variables to the model one at a time* One begins the procedure by postulating a one 
parameter model 

? = 6 , . I 

where is singly the mean of the data. A tableau of candidate independent model 

variables « X2» Xp is available. The tableau for the airplane longitudinal 

equations of motions is given by table 1. 

Next the correlation coefficients 



are calculated for each x- Then the x. corresponding to the 

greatest ^xjy selected to enter the regression* N«>w the model is 


y = 6^ + e.x 


X. 4- £ 


Next, the correlation coefficient for each remaining (i = 2, «•*, j^l# 

j+1, ***, P) is calculated on and y where y = corre- 

lation coefficient is 
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r is read the partial correlation of y on Xj^ given that is in the 

regression* 

The giving the largest value of r ^ is selected and entered Into the 

model a *The nfK>del is now 


y « 8- 


9.x* 
3 1 




Again a partial correlation coefficient, calculated and the 

giving the largest enters the model- 'idiis process continues until the 

remaining candidate variables offer no statistical improvement in the model at a 
level of significance selected by the user* 


At each step of the regression, one has several pieces of information avail-* 
able- First, the most obvious is the model structure* Second, there is a least- 
squares estimate of the parameters for each itodel selected* Ihird, the F-value of 
significance can be calculated by 


F 



S 


Fourth, a partial F value Fp given by 


y _ N -..J ^ 1 

P m - 1 1 - r 

yx, ♦x . 

^ k 2 


is calculated giving the relative statistical significance of each variable in each 
model, given that the othex variables are present* 

The stepwise regression used for this analysis differs slightly from that in 
reference 4. Hence, it is required that the linear model (that model containing only 
a, q, and 6^) be fit first and the nonlinear terms considered next* 'Ihis con- 
straint allows the user to view the linear model performance before searching for a 
more complex model* 

The best model is chosen as the one corresponding to the largest total F-value* 
If that peak is not clear, then the set of is searched on several of the 

better candidate models to choose the most consistent and parsimonious model* 

Simulated data created from a nonlinear model at high angle of attack have borne 
out these methods* With no noise on either the dependent or independent variables, 
the correct model was selected and correct parameter values identified* With noise 
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on all variables (worst case)# the procedure selected the correct model structure in 
the rollings and yawing-moment equations* In the side-force equation, two of the 
nonlinear terms were not selected in the presence of noise; however, no term that did 
not appear in the true model was selected either. Hence the technique selected, at 
worst, a parsimonious model and, at best, the true nonlinear model. 


DATA ANALYSIS 

The data were obtained from a low wing general aviation airplane with a modified 
wing leading edge. Because of the leading-edge modification, the test airplane could 
be trimmed at amgles of attack that were well above the stall angle of the unmodified 
wing. Taking advantage of this large range of trim angle of attack, two major data 
sets were created. 

The first data set consists of the results of 21 transient small amplitude 
longitudinal maneuvers initiated from steady-state flight conditions. The trim 
values for a in these 21 maneuvers varies from 2.3^ to 20. 2®. Ihese small ampli- 
tude maneuvers were analyzed with respect for the normal-force and pitching-moment 
coefficients. The parameters, C 2 ^, and , extracted 

from the 21 small amplitude maneuvers are plotted in figure 1. 

The second data set consists of several maneuvers made from steady-state flight 
trimmed at angles of attack of about 2®. Prom this trim angle, a large elevator 
pulse was initiated, resulting in a large excursion from the trim a. A con^arison 
of such a large control input and its resultant a deviation with a typical small 
control input (from the 21 small amplitude runs) and its a deviation is depicted in 
figure 2. 

The second set of data was then analyzed in two ways. In the first analysis, 

the entire data length from a particular maneuver was analyzed by using the stepwise 

regression. This technique is the same as the one that was applied to the 21 small 
amplitude maneuvers. The results of this analysis for two particular large maneuvers 
(run 18 and run 19) are given in figures 3 and 4. It is seen in the figures that 
trend lines (dashed lines) have been extrapolated from the trim value of that param- 
eter to higher a values. The trend line reflects first derivative values of the 
parameter (such as Cg and Cz 2 ) that were selected by the stepwise regression. 

A con^arison of these trends with the small an^litude maneuver parameter stimates is 

seen in figures 5 and 6. The trend lines linearly interpolate the small ampliti^de 

data very well. Ihe trend information is not chosen for all parameters i*' all 
maneuvers. In such cases, the trend line is represented by a horizontal vrajectory 
signifying a zero slope for that parameter with a. 

To improve upon the trend information or rather to achieve a better resolution 
of the fine structure of the parameter as a function of a trajectories, a second 
method of analysis was applied to the large amplitude data* 

In this second method, the data from each maneuver were divided into bins as a 
function of a. Tliat is, all measured data corresponding to 0 < a < 4* %#ere put 
into the first bin. Drta corresponding to 4®<a<8* were put into the second bin 
and so forth until six bins each of 4* width yrete filled, alien five new bins for 
2® < a < 6* , < Q < 10®, and so forth, were created and filled with the corre- 

sponding measured data. Ohe number of data points and the mean a value of those 
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points are given for each bin for runs 18 and 19 in table 2* Each bin was analyzed 
by using the stepwise regression and the results are given in figures 7 and 8* *lhe 
4^ bin width as chosen as e means of including enough data in each bin but also to 
restrict each bin to as small a range of a as possible. Ihe smallest range of a 
prevents data from several regions of aerodynamic behavior from diluting one another. 
With greater data rates or longer maneuvers it might be possible to make even more 
narrow bins. 

The results of the bin analysis are overlaid with the small amplitude analysis 
and the entire large amplitude analysis in figures 9 and 10. The solid circles 
represent the parauneter value extracted for a given bin and is plotted at the 
mean a for the points in that bin. For certain bins, no linear parameter was 
extracted; hence not all bins have a representative parameter value. However, cer- 
tain bins yielded trend information, as in the analysis of the entire data length. 

The trend lines, where available, are plotted at the mean value of a for the points 
in that bin. One sees that though the entire run analysis can give a general trend 
in the parameter behavior at higher a, the partitioning of the data into bins 
provides a fine structure about that trend. Hence more information can be extracted 
from the Scune data by dividing it into bins than by simply analyzing the entire data 
length. 

The ordinate value at whicn a slope line is plotted in this report is determined 
by extrapolating the ordinate values of previous bins to achieve a reasonable place- 
ment. The tenptation to treat the ordinate value as equal to zero in cases where 
only a slope is identified must be suppressed since zero is not the "best** value for 
the ordinate but rather represents the value roost conpatible with the best overall 
model. Because of correlation among the variables in the adequate model, the ordi- 
nate value need not be included as a pareuneter in the best model in certain cases. 

An example will better explain the previous paragraph. Consider the binned 
results for C„ of run 19 as shown in figure 10(a). If one considers the equation 

q 


16® bin (a = 14.4®) of run 19 would give 


= 0 - 1.82(14.4® - 2.3®) * -22.0 


= 0 + (a - 
Z Z ^ 

q qot 


and, for example, the 12® to 


when such a formula is used, one develops the plot for against a shown in 

figure 11. ^ 

Certainly it is better to use the additional l.iformation provided by the sur- 
rounding bins (as in figs. 9 and 10). However, it is better still to develop an 
analytical expression to incorporate the ordinates and slope information from sur- 
rounding bins into an optimal polynomial fit through all points, obeying the given 
slope information. The development of such an expression is currently being studied. 
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DISCUSSION 


For the normal-force equation, it is seen in figures 9(a) and 10(a) that 
can be extrapolated from -4.3 a . a • 2.3® to -1.3® at a • 20® by sinply using 

2 2 

a C /3a = 2C from the large maneuver entire data length. Ihis general trend 

CL 

corresponds to the values extracted from the small aunplitude transient maneu- 


vers* However, the behavior between these end points is only approximated by s.ich an 
analysis* Ihe finer behavior is washed out by the averaging of all the data* When 
one bins the data, it is seen that the trend from trim to about a = 7® is a 
slow increase, then a rapid increase in C« between a = 7® and 12®, and again a 


slower increase from a = 12® to 20®* A similar behavior is seen for 


More- 


over, for both runs 18 and 19, it is seen that slope information Cr 


qa 


q 

is 


rtracted from certain bins* The slopes are plotted about the mean a for their 
respective bins. In run 18, one sees that the ~>lope crosses the overall trend line 
between a = 7® and 1 1® reinforcing the information from that the system 


behavior is changing most rapidly in this region* 
in slope after its maximum decrease at a 


Ran 19 gives the general decrease 
9.5®* For C 2 g t the entire run 


analysis gives no nonzero Cgg ^ value. This is represented by the zero slope line 

in figure 10(a)* Though the bins of run 19 provide no additional information, run 18 
bins provide a slope Cgg ^ at 8.7® and shows that C 2 g is about constant at -0.8 

from cIq to about a = 8®, and then C 2 ^ stays constant at about -1 from a = 9® 

to 20®. The standard error at high a makes an exact measure of C2r impossible. 

^e 


For the pitching-moment coefficient equation, one should r' * er to figures 9(b) 
and 10(b). The analysis of the entire data length of run 18 failed to give any 
nonzero trend information for any of the pitchingnmoment derivatives* This is noted 
by the horizontal line representing a zero slope trend. The analysis of the entire 
data length for run 19 g trend lines for Cj}|^ and C^* By binning the data, 

one sees the interesting fine structure in run 18 for derivative is 


seen to be reasonably nstant from to about 8®. Between a = 8® and 12®, C* 

drops rapidly, steadies from a = 12® to 15*6® and then rapidly decreases as shown 
by the exiscence of bins between a - 16® and 20®. Though not as 


visible, the run 18 behavior is not contradicted by run 19 binned analysis* 


The variation of 



about the overall trend line for run 19 shows a rapid 


increase in between a == 4® and 10® in that the binned trajectory cross the 

trend line in that region. Ran 16 adds no new inforntation but also does not 
contradict run 19. The derivative is shown, in the binned analysis of run 19, 
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to have a constant value of about -1.95 from Oq to a = 7® after %ihich it 
increases bet%ieen a - 7* and 10* for a new constant value of about this 

behavior « especially in regions of change, is consistent with the behavior of the 
other derivatives in both the binned d*\ta and the 21 small ac^litude maneuvers* 


One final aspect of the application of the stepwise regression for both large 
and small ang>litude maneuvers performed over a large range of trim angles of attack 
is shorn in figure i1. One can take the average excursion in a to be given hy 


N 

i E (&o )2 

N i=1 1 

a given run* 


where is fcr the ith value of in a data string for 

Ihe nvBiber of terms selected for the optimal model xik s^own as a func- 


N 

tion of ^ (Aa^)^ and ^or Uie 21 small amplitude maneuvers and the large 

amplitude maneuvers* 


In the semilog plot in figure 12, points where linear models suffice are repre* 
sented by circles* Maneuvers %#hich required nonlinear models are denoted by squares* 
The number in the circle or square is the number of terms present in the adequate 
model* A line separates the plane into two sections* Below the line one finds that 
a linear model is generally adequate* Above the line, generally a nonlinear model is 
required. From the plot one sees that justification of a linear model depends both 
on the excursion from and the behavior of the system in the neighborhood of a^* 

Hence even at large a^, a linear model suffices for an airplane with a reasonably 
linear lift curve at that a* Also, at low o^, a nonlinear model is required if the 
excursion from Oq is large* 


CONCLUDING REMARKS 

It has been shown that information regarding the behavi';>r of an airplane at high 
angles of attack exists in data that results from large ang>litude longitudinal maneu- 
vers from a low trim angle of attack* A technique has been demonstrated which allows 
for the determination of overall trend lines indicating the mean behavior of the 
longitudinal linear stability and control derivatives over a large range of angle of 
attack. Though this trend line can be useful for gross estimates of airplane behav- 
ior at high angle of attack a, a technique of partitioning the data string was shown 
to give some fine structure and detail about the trend lines* These detailed esti- 
mates from partitioning were shown to be reliable in that they agree with estimates 
using small perturbation from trim condition in certain regions* 

The partitioning is especially useful to detect regions of rapid change in 
parameter values* The regions are usually indicated by the significant contribu. ; on 
of one or more nonlinear terms into the model* These transient regions are areas 
that should be more closely investigated by further flight tests. They could be 
interesting areas of nonlinear behavior involving multiple states or hysteresis* 
Future research should include means of extracting still more information from these 
"transient data" regions* 

The modified stepwise regression, as presented, has been shown to be reliable in 
choosing an adequate model* Because it does not test all possible models, one is not 
assured that he has the absolute best model* However, work with simulated nonlinear 
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■odels has shorn that the model chosen, even under noisy flight conditions is ade- 
quate* Future wrk should include development of an information criterion for a best 
model* 


Finally# it was shown that in some regions the linear model is the best model# 
but that over a broader range of flic^t conditions# a nonlinear model is required* 
Moreover in certain regions of trim angle of attack even for small perturbations 

from trimmed flight# a linear model is inadequate for stability and control analysis* 

*Ihe techniques presented in this paper show promise for the task of model struc- 
ture determination and analysis of flight data from transient flight regimes* Itie 
future development of information criteria and experience developing a better set of 
candidate model variables will enable the user to better analyze data from transient 
maneuvers such as are found in prestall and poststall regimes* 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton. VA 23665 
August 17# 1981 
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APPENDIX 


EQUATIONS OF MOTION 

The following equations of motion are given in this appendix. 


• ^ PV^S 

u » -qu rv - g sin 0 - C 

2m X 


* a pV^S 

V = -ru + pw + g cos 9 sin ^ + ■■ ■ C 

2xn Y 


PV^S 


w = -pv qu + g cos 6 cos ^ 

2m Z 


,^Y ~ 1 . ^XZ, . . PV^Sb „ 

P = qr\— ; I + — <pq + r) + 




X / X 


\ - V, 2 

— ; + T — <r - 


q = pr — ^ ; 

\ -Y / ‘y 


2Iy m 


^X ~ ^Y| ^XZ • . pV^Sb 

r = ?q|-^ 1 + — (p - qr) + 

z / z z 


9 = q cos - r sin 


= p + (q sin 4> + r cos tan 9 


/ 2 2 2 
V = /u ^ V t w 


w 

a = tan — 
u 


• w 

a « — 
u 


6 = sin 


-1 V 


• V 

0 « - 

V 
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where 


C = C + c (a - o ) + C ^ + c_ (6 - 6 ) + c (o - 

Z Z,o 2 o Z 2V Zr e e,o Z o 2V 
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Z 2 o Zf o e e,o Z^2 o 

at o QL p 

e 
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m ^2 c m 3 o m 4 o mb o 

B a a a 


6 7 8 

^ C* ^(a - a ) + C* _(a - a ) + C* _(a - a ) 

m 6 o m7 o m 8 o 

a a a 


where 


♦ C 


C* “ C + C ( T — C + cos 0 ] 
m,o m,o m«l 4m Z,o _„2 o 

“\ 2v y 


c = C + 7 -^ c c 
m m 4m m* Z 
a a a a 


C = C + C [1 + 7 ^ C 
m m m* \ 4m Z 
q q a\ q 


pSc 

r* = c + C c 
m. mj, 4m m» Z. 
60 a 0 

e e e 


pSc 

C' = c + c ^ ^ 

m m m* 4m Z 

qa qa a qa 
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pSc 

c' = c + c c 

F 2 m2 m* Z 2 

a a a a 


C’ = C + C ^ — c 

m . m* 4m 

^ a 0 a a C a 

e e e 


pSc 

^ ^ = C ^ + C ^ C„ ^ 

mo2 mrt2 m* 4m Z^2 

3 a 6 a a 6 a 


PSc ^ 

C = P 4- P ^ C 

*b 2 nip2 m* 4m Z^2 


C . = C . + C C . 

mi mi m* 4m Z i 
a a a a 


Assuming steady-state \V^ = Pq === qo ~ ~ initial flight conditions the 

longituc Inal equations become 


-qw - g sin 8 + 


n 

qu + g cos 9 + ~ - C 

2m Z 


• pV^ Sc _ . 


and the longit' iinal output can be written as 


= “(u + qw + g sin 9) 
X g 


a„ = - qu - g cos 6) 

Z g 
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and for the equation error form, the aerodynamic coefficients can be written as 


2gm 

pV^S 


C 


X 


2gm 

PV^S 


a 


z 


c 


z 


21 . 


pv^sc 


C 

m 
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Figure 1.- Aerodynamic coefficients 
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(b) t>itching-moment coefficients. 
Figure 1.- Concludede 
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(a) 2-force coefficients. 

Figure 3.- Aerodynamic coefficients from run 18 (large amplitude)^ 
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(b) Pitching-moment coefficient. 
Figure 3.- Concluded. 
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(a) Z-foice coefficients. 

Figure 4.- Aerodynamic coefficients from run 19 (large amplitude). 
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(b) Pitching-noment coefficients. 

Figure 4.- C' '"eluded. 
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O Small amplitude maneuver 




(a) Z-force coefficients. 

Figure 5.- Comparsion of aerodynamic coefficients from run 18 and small 

an^litude maneuvers. 
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(b) Pitching-moment coefficients. 
Figure 5.- Concluded. 
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(a) Z-force coefficients. 

Figure 7.- Aerodyneunic coefficients from binned run 18. 
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(b) Pitchinq^moroent coefficients. 
Figure 7.- Concluded. 
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(a) Z-force coefficients. 

Figure 8.- Aerodynamic coefficients from binned run 19. 
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(a) Z-force coefficients. 

Figure 9.- Comparison of aerodynamic coefficients from run 18 and small 

amplitude maneuvers. 
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Figure 10*- Comparison of aerodynami 
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(b) Pitching-moment coefficients. 
Figure 10.- Concluded. 
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Figure 11.- Run 19 binned slope values for Cg referenced to trim a. 
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